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TITLE OF THE INVENTION 



DETERMINATION OF A SELECTED POSITION OF A MOBILE COMMUNICATIONS DEVICE 

IN A COMMUNICATIONS NETWORK 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] This application is based on and hereby claims priority to German Patent Application 
No. 103 18 590.9 filed on April 24, 2003, the contents of which are hereby incorporated by 
reference. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

[0002] The invention relates to determining a position of a mobile communications device in 
a communications network (localization). 

2. Description of the Related Art 

[0003] With the increasing spread of mobile communications, the demand for additional 
services in the mobile radio system is also growing. 

[0004] Location Based Services in this case are taken to mean additional services of mobile 
radio operators which are offered or provided to users of the mobile radio system on a location 
basis, i.e. depending on a position or a geographical location of relevant users, for example 
location or distance-dependent usage tariffs or orientation aids for search and rescue services. 

[0005] Fundamental to a "Location Based Service" is therefore the localization or position 
determination of the user concerned or of their mobile communications device. 

[0006] Different techniques are known for such localization of mobile communications 
devices in communications networks, for example determining the position on the basis of delay 
time determination or delay time measurements of communications signals of a mobile wireless 
communications device to a base station of a communications network (Rappaport T.S.,Reed 
J.H. et al., "Position Location Using Wireless Communications on Highways of the Future", 
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IEEE Communication Magazine, pp. 33-41, Oct. 1996, DE 198 36 778 Al) or a localization by a 
satellite-assisted system such as GPS. 

[0007] The delay time-based position determination known from DE 198 36 778 Al is 
executed for a mobile telephone, generally a mobile station, in a GSM communications network 
(= Global System for Mobile Communications) (Eberspacher, J.; Vogel, H.-J.: GSM. Global 
System for Mobile Communication. Stuttgart, Leipzig: Teubner, 1999; Jung. P.: Analyse und 
Entwurf digitaler Mobilfunksysteme. Stuttgart, Leipzig:Teubner, 1997; Kennemann, O.: 
Lokalisierung von Mobilstationen anhand ihrer Funkmessdaten. Nummer 11 in Aachener 
Beitrage zur Mobil- und Telekommunikation. Aachen: Verlag der Augustinus Buchhandlung, 
1997) in accordance with a TDMA (Time Division Multiple Access) mobile radio technology. 

[0008] An individual mobile station, which has registered with a fixed local base station (a 
base station which is controlling the call) is currently assigned a free time slot in a TDMA frame 
by the base station. 

[0009] The communications signals intended for the mobile station concerned go to this time 
slot in signal packets, known as bursts, with a length of 15/26ms, from the base station or the 
communications signals or bursts sent from the mobile station must arrive at the base station. 

[0010] As a result of scattering the communications signals emitted by the base station find 
their way to the mobile station via different paths (multipath propagation), in which case they are 
attenuated depending on their frequency. 

[001 1] A receive field strength of the communications signals received by the mobile station 
is thus not only dependent on the distance from the mobile station to the base station, but also 
on the frequency and topographical circumstances between mobile station and base station. 
Therefore the individual data packets are sent on different carrier frequencies, which allows 
selective interference on one frequency to be distributed between a number of users. 

[0012] However a precise synchronization between mobile station and base station is 
needed to enable this to be done. This synchronization is additionally made more difficult by a 
mobility of the user, since the mobile station is now located at changing distances from the base 
station and its communications signals exhibit different delay times. 
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[0013] To equalize the different delay times and enable frame-synchronized data to be 
delivered to the base station, the mobile station measures the signal delay time to the base 
station and uses this to correct the start of sending its bursts. 



[0014] The signal delay time is encoded in a "timing advance" (TA) and exhibits a 
dependency on the distance between mobile station and base station controlling the call. 

[0015] 64 stages are available for the TA which are (bit) encoded with the values 0 to 63 and 
represent the delay time. 

[0016] Since positions of base stations are known the position of the mobile station can be 
concluded from a TA or from the signal delay time. 

[0017] A measurement accuracy for determining the delay time relates to a bit duration, 
which is 48/13 ps in GSM corresponding to a single path length of approx. 554m. 

[0018] Determining the position of a mobile communications device in a UMTS (= Universal 
Mobile Telecommunication System) network is known from TS 25.305 V3.1.0: stage 2 
"Functional Specification of Location Services in UTRAN" (release 99), 3GPP TSG-RAN-WG2, 
2000. 

[0019] In the corresponding UMTS mobile radio standard on which the UMTS network is 
based, determining the position of a mobile radio device is already explicitly included in the 
standard or is required by it (TS 25.305 V3.1.0: stage 2 "Functional Specification of Location 
Services in UTRAN" (release 99), 3GPP TSG-RAN-WG2, 2000). 

[0020] Further methods for localization of a mobile communications device in a 
communications network are known from United States Patent Number 5,883,598; United 
States Patent Number 6,094,168; United States Patent Number 6,108,553. 

[0021] A non-linear filter, a Prior Density Splitting Mixture Estimator (PDSME) is known from 
Rauh Andreas, "Nonlinear Measurement Update And Prediction: PRIOR DENSITY SPLITTING 
MIXTURE ESTIMATOR", submitted to IEEE Transactions on Automatic Control, December 
2002. 

[0022] This non-linear filter, this PDSME, is based on new Gaussian mix filter algorithms for 
non-linear system state equations and/or non-linear measurement equations. 
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[0023] Usually linearization errors occur if non-linear system state equations and/or non- 
linear measurement equations are replaced by linearizations during the use of linear filter 
techniques. 

[0024] The PDSME reduces these linearization errors by splitting the a-priori densities into a 
number of Gaussian mix components each with smaller covariances. 

[0025] The PDSME can be applied both what is known as a "Prediction Step" for the non- 
linear systems state equations and also for what is known as a "filter step" for the non-linear 
measurement equations. 

[0026] Furthermore a measurement is known for the linearization error from Rauh Andreas 
which is comparable with a Kullback-Leibler distance. In addition a desired estimation quality 
can be set for the PDSME. By limiting the maximum number all mix components for the 
PDMSE an upper limit can be specified for the computing effort. 

[0027] The localization methods mentioned have the disadvantages that the positions of the 
mobile communications devices determined by them are inexact and consequently affected by 
large uncertainties. More precise methods however require expensive additional devices and 
costly modifications to the communications network or networks and communications devices. 

SUMMARY OF THE INVENTION 

[0028] An object of the invention is thus to make possible a localization of a mobile 
communications device in a communications network which is accurate and is affected by the 
fewest possible uncertainties, which can be implemented as simply as possible and at low cost. 

[0029] With the method for determining a selected position of a mobile communications 
device in a communications network with at least one first base station set up for a first 
communication with the mobile communications device by first communications signals, a non- 
linear communications model with a deterministic component and stochastic component is 
used. 

[0030] The non-linear communications model includes a deterministic and a stochastic 
component, 

- the deterministic component of which describes a dependency between communications 
signals of a base station and positions of the mobile communications device and 
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- the stochastic component of which describes the uncertainty of the deterministic 
component. 

[0031] By using the first communications signal belonging to the selected position of the first 
base station and the non-linear communications model, the selected position is determined, 
with a non-linear Bayesian filter technique being applied to determine the position. 

[0032] By applying the non-linear Bayesian filter technique to the communications model or 
to the communications signal, a possible location area of the selected position is determined 
from which the selected position is finally determined. 

[0033] The non-linear Bayesian filter technique in the invention is generally taken to mean 
the technique described below: 

- The Bayesian filter technique corresponds to a continuation or modification over time (also 
computing the change over time) of known probability distributions of the system states and 
system outputs by the presence of new measured values. 

[0034] The inventive localization method is based on the idea of obtaining from 
communications signals between at least one base station (transmitter) and one mobile station 
(receiver) available in communications networks such as for example WLAN, GSM, DECT 
networks, distance-relevant parameters and from these geographical information, in this case a 
possible location or distance area or area of uncertainty of the mobile station. 

[0035] It should be noted that technologically the transmitter and receiver of the 
communications signals can also be implemented the other way round. The inventive mode of 
operation is not affected by this. In this case the base station would be the receiver and the 
mobile station the transmitter of the communications signals discussed. 

[0036] The inventive localization method is a further based on the knowledge that the 
communications signals received by the mobile station are dependent on the geographical 
position of the mobile station in the communications network or in relation to the base station. 

[0037] This dependency between at the received communications signals or between a 
signal strength of the communications signals received (by the base station) and the 
geographical position of the mobile station can be described by what are referred to as 
measurement equations, i.e. by the communications model. 
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[0038] These measurement equations or the communications model include a deterministic 
component describing a dependency between the received communications signal and the 
position of the mobile station, and a stochastic component describing an uncertainty of the 
deterministic component. 

[0039] This uncertainty can for example relate to the communications signal and/or to the 
dependency mentioned (model uncertainty). 

[0040] The non-linear Bayesian filter technique applied in accordance with the invention to 
these measurement equations or the communications model for determining the possible 
location area of the mobile station or the selected position significantly improves the quality of 
the localization compared to usual filter techniques such as an Extended Kalmann Filter (EKF). 

[0041] The non-linear Bayesian filter technique can be undertaken once or iteratively a 
number of times in a "filter step". 

[0042] In this possible location area or area of uncertainty, described using stochastic 
variables, such as a mean and/or a variance, the mobile station or its position is finally 
presumed. 

[0043] The presumption can be formulated mathematically by a characteristic value of the 
possible location area or area of uncertainty such as a focal point or an expected value. This 
can then be used as an estimation for the selected position of the mobile station. 

[0044] A particular advantage of the invention lies in the fact that the localization is conducted 
on the basis of communications signals and known positions of base stations which occur in a 
normal operation of a mobile radio system and are available there. This enables expensive 
modifications and expansions as well as additional measurements of existing mobile radio 
systems or for existing mobile radio systems to be dispensed with. 

[0045] The developments described below relate to both the method and to the configuration. 

[0046] The invention and the developments described below can be implemented both in 
software and also in hardware, for example by using a specific electrical circuit. 
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[0047] Further the realization of the invention or of a development described below is 
possible through a computer-readable storage medium on which a computer program is stored 
which executes according to the invention or development. 

[0048] For communication in a communications network between a mobile communications 
device (mobile station), for example a mobile telephone and a base station, for example a 
circular antenna or a circular radiator or one or more sectoral antennas, data, the 
communications signals, is transmitted in signal packets known as bursts. 

[0049] Based on or using at the transmitted communications signals or signal packets 
various parameters relevant to distance can be determined which in their turn can be included 
as a basis for determining the possible location areas or distance areas. 

[0050] Such a distance-relevant, i.e. distance-dependent parameter is for example a field 
strength of a signal packet. 

[0051] The field strength exhibits a natural dependency on the distance between the mobile 
station and the base station (controlling a call) and consequently delivers information about the 
possible location area or distance area (area of uncertainty) of the mobile station. 

[0052] This dependency between field strength and distance can be described by physical 
models which describe a propagation behavior of signals. 

[0053] For the determination of the selected position a characteristic value of the possible 
location area or area of uncertainty can be determined, such as a focal point or an expected 
value which is then used as an estimation for the selected position of the mobile station. 

[0054] The quality of the localization can be further improved by using a user model which 
describes a movement of the mobile communications device. Thus such a user model can limit 
the maximum step length or movement distance for a prespecified time step. 

[0055] The invention is especially suitable for use in at the environment of a digital cellular 
mobile radio system such as a GSM network, and there for example for localization of a GSM 
telephone (mobile telephone). 
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[0056] In this case when the invention is used only the data available to a mobile telephone 
will be used, in which case costly changes to not have to be made either to the GSM network or 
to the mobile stations in the GSM network. 

[0057] For example the positions of the individual base stations and their antennas as well as 
their characteristics which provide information about the coverage area of the relevant antenna 
are known from a GSM network. 

[0058] The mobile telephone for its part always stays in contact with the antennas that can be 
received for a correct connection set-up so that it can be assigned by the network the antenna 
which is most suited for a call. To this end it measures the field strengths of the receivable 
antennas as well as specific signal delay times which are then also known. 

[0059] On the basis of this available information the mobile telephone is then localized in 
accordance with the inventive method of operation. 

[0060] The invention is also suitable for use in the environment of other digital cellular mobile 
radio systems such as a WLAN or of a DECT network (11], and there for example for 
localization of a DECT mobile telephone. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0061] These and other objects and advantages of the present invention will become more 
apparent and more readily appreciated from the following description of an exemplary 
embodiment, taken in conjunction with the accompanying drawings. The exemplary 
embodiment is in this case subdivided into a basic part which presents the fundamentals of the 
inventive method of operation and a related part which clarifies or specifies on the basis of a 
concrete numerical example the inventive method of operation as well as the results. In the 
Figures, 

Figure 1 is a graph in which it an assumed linear approximation for the non-linear 
communications model is presented; 

Figure 2 is a graph in which the uncertainty of the communications signal through noise 
is presented; 

Figure 3 is a block diagram providing an overview of the PDSME algorithm with a 
linearized filter step (upper section) and a prediction step (lower section); 
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Figures 4A and 4B are graphical representations of an evaluation scenario for 
localization in accordance with the exemplary embodiment; 

Figures 5A and 5B are graphs in which a localization using an extended Kalman filter is 
shown (first filter step Fig. 5a; 170th filter step Fig the. 5b); 

Figures 6A and 6B are graphs in which a localization using the PDSME in accordance 
with the exemplary embodiment is shown (first filter step Fig. 6a; 170th filter step, Fig. 6b; 

Figures 7A and 7B are graphs in which results of the EKF filtering and the PDSME 
filtering in accordance with the exemplary embodiment in relation to the correct results are 
shown; 

Figure 8 is a graph in which a number of Gaussian mix components which depend on 
filter steps are shown. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

[0062] Reference will now be made in detail to the preferred embodiments of the present 
invention, examples of which are illustrated in the accompanying drawings, wherein like 
reference numerals refer to like elements throughout. 

Exemplary embodiment Localization of a DECT mobile telephone in a DECT network 
based on a non-linear filter technique 

Basics 

[0063] The application of a non-linear Bayesian filter technique in the localization of a mobile 
radio telephone is described below. The application is described using as an example a 
localization of a DECT mobile telephone in a DECT network with a number of base stations. 

[0064] Signal strengths (as field strengths of the communications signals received by a 
mobile telephone and measured there (of a base station transmitting the signals) are dependent 
on the position of the mobile telephone in relation to the base station sending out the 
communications signals. This relationship is described through non-linear measurement 
equations (non-linear communications model). 

[0065] These non-linear measurement equations or this non-linear communications model 
includes a deterministic component which describes the signal strengths received as a function 
of the position and a stochastic component which takes account of model errors and 
measurement noise. 
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[0066] In addition user models are also considered which bring with them knowledge about a 
(spatial) movement of users of the mobile telephones. 

[0067] A new non-linear filter technique, known as a "Prior Density Splitting Mixture 
Estimator" (PDSME), assignable as a Gaussian mix filter algorithm significantly improves a 
localization quality compared to standard filter techniques such as an "Extended Kalman Filter" 
(EKF) (cf. Figures 5 to 7). 

[0068] Typical applications for localization tasks are determining the positions of mobile 
communications devices in WLANs, GSM networks or DECT networks. 

[0069] For position determination information should be able to be used in such cases which 
is already available during of a normal operation phase, i.e. for "normal" communication 
between mobile device and base station. 

[0070] The localization method known from T. Roos, P. Myllymaki, and H. Tirri, "A Statistical 
Modeling Approach to Location Estimation", IEEE Transactions on Mobile Computing 1, pp. 59- 
69, January-March 2002 for mobile telephones in GSM networks is based on location- 
dependent signal strengths of the communications signals sent out by base stations, able to be 
received and measured by the mobile telephones. Because of the location dependency these 
are characteristic for the position of the signal recipient, in this case for the position of the 
mobile telephone or of a user of the mobile telephone 

[0071] During the normal operation phase the received signal strengths of all receivable 
communications signals of the relevant transmitter (base stations) are measured and compared 
by the mobile telephone. This is how a handover of the mobile telephone between different 
transmitters or base stations is controlled. 

[0072] Based on these measurements made during normal operation localization methods 
can be developed in order to additionally determine the position of the mobile telephone in 
relation to the transmitter or in the GSM network. 

[0073] Further procedures are known from A. M. Ladd, K. E. Bekris, G Marceau, A. Rudys, 
D. S. Wallach, and L. E. Kavraki, "Using Wireless Ethernet for Localization", in Proceedings of 
the 2002 IEEE/RSN Intl. Conference on Intelligent Robots and Systems, pp. 402-408, EPFL, 
(Lausanne, Switzerland), October 2002 and S. Panzieri, F. Pascucci, and G. Ulivi, "An Outdoor 
Navigation system Using GPS and Inertial Platform", IEEE/ASME Transaction on Mechatronics 
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7, pp. 134-142, June 2002 for localization of mobile communications devices. These 
procedures for localization are based on an application of non-linear status estimations. 

[0074] The application of a new non-linear filter technique for the localization of DECT mobile 
telephones will be described below. 

[0075] To make it possible to localize mobile telephones, as is known from B. H. Fleury and 
P. E. Leuthol, "Radiowave Propagation in Mobile Communications: An Overview of European 
Research," IEEE Communication Magazine 23(2), pp. 70-81, 1996, a stochastic approximation 
model for the propagation of radio waves or receivable signal strengths of communications 
signals in a localization environment is identified or developed. 

[0076] A corresponding propagation model can also basically be obtained through physical 
propagation models of electromagnetic waves. 

[0077] An adaptation of purely physical propagation models to real (localization)environments 
is however very complicated since reflections or interference or other physical phenomena 
influence and falsify the assumed physical propagation. 

[0078] Since real propagation parameters which describes the electromagnetic (propagation) 
characteristics in the real localization environment in which the mobile telephone is to be 
localized are only partly known or determinable the use of such physically-based propagation 
models for such real applications is excluded. 

[0079] Therefore a propagation model which is based on actual measurements of the 
(logarithmized) received signal strengths of the communications signals in the actual localization 
environment or which can be derived from it is to be preferred over theoretical physical 
propagation models. 

[0080] The development of such a measurement model or communications or propagation 
model derived from actual measurements requires a calibration within the framework of which 
parameters of the measurement model are defined. 

Calibration ("model generation phase") 

[0081] For such a calibration the receiver (used in the sense of a measurement device) i.e. 
the mobile telephone, is moved to a number of positions within the localization environment, for 
example at grid points obtainable by placing a grid over the localization environment. 
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[0082] At each grid point the field strengths of all base stations which can be received at this 
point are measured and stored in a map of the localization environment. In this way the field 
strength distribution of each base station over the localization environment is determined. The 
field strength distributions can be modeled by the measurement model. 

[0083] Based on this model information (and on field strengths measured at an actual 
position) the actual position assumed by a receiver or by the mobile telephone in the localization 
environment can then be determined (localization in a localization phase). 

[0084] Furthermore a description is given below of how a stochastic measurement model can 
be identified in the model generation phase. 

[0085] This stochastic measurement model is an analytical, deterministic measurement 
function which describes the logarithmized received signal field strength as a function of position 
co-ordinates. In addition a stochastic component is identified which takes account of model 
uncertainty and also measurement noise. 

[0086] Such a stochastic approach to the localization takes account of the uncertainties in 
the position determined through probabilities or probability density functions. 

[0087] As well as the measurement model for the receivable communications signals the 
quality of the localization is further improved by the fact that a user model is included as well. 
Via this user model knowledge about a (spatial and/or temporal) movement of a user can be 
incorporated. 

[0088] Through such user models maximum "step lengths" of users in a time step which is 
defined via two associated or corresponding signal measurements (at the start and at the end of 
the times step) can be restricted. 

[0089] An additional improvement of the localization quality can be achieved by recursive 
position determination through a combination of multiple different consecutive measurements. 

[0090] The new stochastic approach used here for localization results in non-linear multi- 
dimensional measurement equations. An exact solution of the Bayesian filter problem needs to 
complicated non -Gaussian probability densities which describe the position sought along with 
its uncertainties. 
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[0091] These probability density (functions) can be approximated by Gaussian mix densities 
K. Ito and K. Xiong, "Gaussian filter for Nonlinear filter Problems/ 1 IEEE Transactions on 
Automatic Control 45, pp. 910-927, May 2000; D. Alspace and H. Sorenson, "Nonlinear 
Bayesian estimation using Gaussian sum approximations," IEEE Transactions on Automatic 
Control 17, pp. 439-448, 1972. 

[0092] This approximation is undertaken here by a newly developed "Prior Density Splitting 
Mixture Estimator" (PDSME) which has been developed to make it possible for the user to set 
the quality of the localization - by contrast with usual, non-linear, known filter techniques such as 
an Extended Kalman filter (EKF). 

[0093] However specifying an upper limit for a computation effort by restricting the number of 
Gaussian mix components cannot reduce the localization quality. In addition, by the 
representation of densities in the Gaussian mixtures a recursive method of operation both for 
the non-linear predictions steps and also for a non-linear filter steps is possible. 

[0094] Furthermore the first part describes how the stochastic localization problem can be 
formulated. The deterministic components and the stochastic components of the measurement 
model will be explained below. The use of Bayesian filter algorithms based on a Gaussian mix 
approximation for the exact probability densities used in the new PDSME will be described. 

[0095] In the second part results of a localization of a DECT mobile telephone in a DECT 
network will be presented and the quality of the localization using the PDMSE will be compared 
to that using a known EKF. 

Formulating the localization problem 

[0096] The localization task for a mobile communications device in a communications 
network can be divided up into two basic problems: 

1) the identification of the measurement model with a deterministic and a stochastic 
component in the model generation phase; 

2) the determination of the probability density of the position sought in a Bayesian filter step 
in the localization phase. 
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[0097] The measurement model 
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of the #h transmitter, ^ = 1, ...,N, as a deterministic non-linear function h ^ (xk) of the position 
coordinates xk=[xl,k, x2,k] for each of the N different transmitters. 

[0098] The stochastic component vk of the measurement model is represented by an 
additive uncertainty. Additive uncertainties in the measurement model of the logarithmized 
receive signal strengths correspond to multiplicative uncertainties of the noise signal power P 
which are to be attributed to an influence of unspecified attenuations. 

[0099] Each transmitter can be identified by its own ID which is transmitted by the base 
station during a communication. 

[0100] This makes the localization problem considerably simpler since each measured value 
can be directly assigned to the appropriate non linear measurement equation. 

[0101] In the localization phase the non-linear filter algorithms will be applied in each time step 
k. However without the additive uncertainties vk the implicit, non-linear relationships between 
the vector with the measured logarithmized received signal strengths yk of the transmitter can 
only be resolved numerically for example by a recursive least squares method. If the 
uncertainties are also taken into account the measurement equations can be used to update the 
position sought in the filter step. 

[0102] The further described method of operation for the localization is based on an 
approximation of the exact solution of the Bayesian filter step by Gaussian mixture densities 
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[0103] with L components, defined by non-negative weights or expected values p and 
covariances C. 

[0104] If the measured value y is available the exact solution of the Bayesian filter step is 
produced: 

f x e (x k |Yk, ;) = c f x p (x k ) fv, c(Yk, C - hC(x k )) 

[0105] with the non-linear measurement function h(x), the additive uncertainty density 
functions f (v) and a normalization constant C. 

[0106] A-priori knowledge about the position sought is represented by the density function f(x). 

[0107] For each new measurement the previous a-posteriori density f(x|y) is interpreted as a 
new a-priori density so that the position sought or determined can be recursively updated. 

[0108] Similarly to the Bayesian filter step a prediction step which is used to describe the non- 
linear user model 

*k + i = § k fek) + w K (2) 
is likewise executed for approximation of the exact densities by Gaussian mixtures. 

[0109] The non-linear function a(x) is a deterministic model of a user movement. 
Uncertainties are again taken into account by additive noise. 

Stochastic modeling 

[0110] The measurement model is identified using measurements of the logarithmized signal 
strengths for each transmitter at grid points of a grid which covers the localization environment. 
The deterministic (h(x))and the stochastic component (v) of the measurement model (1) are 
identified before the localization. 

[01 11] The deterministic, analytical component h(x) is a measurement equation which 
describes the received signal power as a function of the position coordinates. The stochastic 
component v is a model for the uncertainty of the deterministic component. These uncertainties 
are made up of a spatial uncertainties corresponding to model errors and measurement noise 
over time. 
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Deterministic measurement model of the logarithmized received signal strengths 

[0112] Fig. 1 shows that the acceptance of the logarithmized receive signal strengths over a 
distance of a number of meters is almost linear. In other documents similar assumptions in 
relation to the measurement model are referred to as "linear-loss-model" J. Lahteenmaki, 
"Radiowave Propagation in Office Buildings and Underground Halls", in Proc. 22nd European 
Microwave Conference EurMC'92, pp. 377-382 (Espoo, Finland), 1992.or "linear-scope-model 
B. H. Fleury and P. E. Leuthol, "Radiowave Propagation in Mobile Communications: An 
Overview of European Research," IEEE Communication Magazine 23(2), pp. 70-81, 1996. 

[01 1 3] Mathematically linear acceptance of the logarithmized receive signal strengths for two- 
dimensional position coordinates x can be described by N independent measurement 
equations. 



[0114] The parameters m£ and P£ of the positive semi-definite quadratic form and the additive 
offset A£ are identified for each transmitter 1 N. 

[0115] The 6N parameters m£, P£ and £C are determined by a method of recursive least 
squares in which the deviation between the approximated measurement model and the 
measured logarithmized signal strengths at the grid points measured during the model 
generation phase are minimized. To reduce the measurement noise over time the mean value 
over a number of measurements at each grid point Is used to calculate the average of the 
measured values. 

Stochastic modeling of the uncertainties 

[01 16] The stochastic uncertainty model takes account of both the deviation between the 
approximated model described above and the true distribution of the logarithmized receive 
signal strengths measured at the grid points and the measurement noise over time. 

Model uncertainties 

[0117] In this localization approach the deviation between true logarithmic receive signal 
power and the deterministic measurement model h(x k ) is described for each transmitter by a 

Gaussian noise density with a mean value of ju™ and standard deviation cr^ . This Gaussian 




(3) 
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uncertainty is a representation of the approximation error of the deterministic components of the 
measurement model over the entire localization environment. It first represents the mean 
deviation over the localization environment between the model and the true receive power on 
the basis of the incorrect assumptions for the deterministic part of the measurement equation. 
The uncertainties produced from large deviations between the true receive signal power and the 
measurement model are shown in Fig. 1. Secondly measurements have also shown that local 
deviations exist which could have been caused by reflections, non-homogeneous propagation of 
the radio waves and interference. These spatial variations of the logarithmic receive signal 
power can be seen from Fig. 2 in which measurements have been collected at a spacing of 2 
cm. 

Measurement noise 

[0118] In addition to spatial variations of the receipt signal power measurement noise over 
time can be determined by evaluating a series of several different measurements at a fixed 
point. This measurement noise over time is also shown in Fig. 2 for three different measure- 
ments of the logarithmic receive signal power at each measurement point. This noise is in its 

turn approximated by a Gaussian density with a mean of ju™ and standard deviation cr™ . 

Combination of modeling certainties and measurement noise 

[0119] To obtain a simple model with the two uncertainties described above it is further 
assumed that the two uncertainties are independent. They can thus be modeled by a single 
Gaussian density f^ g (u g ) for each transmitter 

1,... t N, which is defined by the mean 
and the standard deviation 

[0120] This model implies not only independent uncertainties but it is also assumed that the 
two uncertainties can be described without taking into account any positional dependency. 
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Filter algorithm 



[0121] This section gives a brief overview of the PDSME (Prior Density Splitting Mixture 
Estimator) used for the localization of radio communications devices. In addition an adaptation 
of the measurement equations derived in Section 3 h(X|<) is introduced to simplify the 
computation or the PDSME estimator. Furthermore a prediction step for a simple user model is 
presented. 

The PDSME measurement updating 

[0122] The PDSME algorithm for the localization of radio communications devices shown in 
this step is based on the calculation of a linearized measurement update for Gaussian mixture 
densities. The measurement updating step of this filter algorithm is shown in a block diagram in 
the upper part of Fig. 3. Splitting is based on the calculation of the linearization error 

for each component of the Gaussian a-posteriori mixture density. This criterion is very similar to 
the Kullback-Leibler distance S. Kullback and R. A. Leibler, "On information and Sufficiency", 
Annals of Mathematical Statistics 22, pp. 79-86, 1951 



between the exact a-posteriori density f * and its approximation(equation) by 

replacing the non-linear measurement * fa f h(X k ) by its linearization ~v3ifc/ a t 
the mean of the ith component of the a-priori-density function. By calculating the linearization 
error (4) the Gaussian a-priori mixture components which contribute most to the approximation 
error of the a-posteriori density are identified. To reduce this linearization error these Gaussian 
mixture components are replaced by being split up into a number of mixture components with 
smaller covariances using splitting libraries which have been optimized offline (see Fig. 4). 
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[01 23] After this "analytical resampling" of the a-priori densities, the filter step can be 
calculated by a bank of EKF, linearization of the measurement equation at the mean of each 
component of the representation of the Gaussian mixture of the a-priori density. 

[0124] Thereafter redundancy is reduced by a merging step by combining a number of 
Gaussian mixture components into a single Gaussian with negligible approximation error. 

[01 25] Adaptation of the measurement model for simplifying the linearization error criterion 

[0126] For Gaussian measurement noise the linearization error (4) can be calculated 

analytically as a linear combination of moments of the densities f* (2fcl#fei<) f or polynomial 
measurement equations h£(X k ). In this subsection an adaptation of measurement equation (3) 
is thus derived. After a few algebraic transformations the quadratic measurement equation 

-*M ' ^ ' (5) 

[0127] can be rewritten with a polynomial function with a modified "measured value" Z ki < and 
the transformed uncertainty U<. Because of the non-linear transformation of the random 

variable u< the probability density function ^ Vf ^ is no longer Gaussian. In the localization 
experiment in Section 5 the exact moments of first and second order 0^ will be calculated to 

determine a Gaussian approximation of ^ v% ^ . It should be noted that the moments U c 
depend on the measured value y k ^. They must thus be recalculated for each new measurement 
of the logarithmic receive signal power and are not unchangeable over time like the parameters 
described in Section 3 p Ui? and a u ,<. 

User modeling by PDSME prediction step 

[0128] Like the filter step the PDSME can also be applied to the calculation of non-linear 
prediction steps. In a similar way to non-linear filter steps the calculation of an approximated 
prediction step also includes an evaluation of a linearization error, a bank of linearized prediction 
steps and the reduction of the number of Gaussian mixture components in a merging step (see 
lower part of Fig. 3). In this work only a linear user model is taken into consideration. The 
predictions step can thus be calculated analytically since a-posteriori density has been 
approximated by Gaussian mixture density in the filter step. 
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[0129] The prediction model uses the linear status equation 



with the mean p w of the additive system noise u>k representing knowledge about possible 
directions and average step lengths of the movement of the user. The covariance matrix C w of 
u)k specifies an estimation for the distribution of the step lengths of the user. For each 
component i = 1.....L of the Gaussian mixture density the predicted Gaussian mixture 
component is then described by the mean 

and the covariance 

The estimated position can then be calculated as weighted overlaying 




of the mean values of all Gaussian mixture components. 
Localization experiment 

[0130] In this section a real localization experiment for the validation of the localization 
approach described is presented for DECT mobile telephones. In these subsections the 
superior performance of the PDSME by comparison with the EKF (Extended Kalman Filter), a 
widely-used standard approach for estimating the status of non-linear systems, is shown. 

Evaluation scenario 

[0131] In this localization experiment N = 10 transmitters have been placed in an internal 
space of approximately 30 m x 30 m on a floor of a building. In the model generation phase the 
logarithmic receive signal part of each transmitter has been measured on a grid with 1 m 
spacing between grid points. IN = 10 measurement equations fy(x k ), consisting in total of 60 
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parameters for the deterministic components and 20 parameters for the stochastic components 
ilk have been identified. Apart from just a few areas which were influenced by the high 
attenuation of the radio waves as a result of steel and concrete walls, the assumed model is an 
appropriate approximation of the distribution of the receive signal power. 

[0132] The initial probability density of the position is selected as a Gaussian density with the 
initial mean value 

«-B£3 

and initial covariance 

02=8 [T is*] mS 

f 

i.e. there is almost no prior knowledge available about the position. 

[0133] In the localization phase a measurement of each receivable transmitter is collected 
along a line, as shown in Fig. 4, with a spacing of Ar = 0,5 m between the true measurement 
points Pi, P 2 , Pi 7 . The position co-ordinates of the measurements have been determined to 
compare the basic actuality with the estimated positions. For the localization of the receiver the 
exact measurement positions have not been used either in the localization approach with the 
PDSME or with the EKF After the measurement update for each receivable transmitter at a 

fixed point Phi = 1 17 has been calculated, a prediction step in accordance with Subsection 

4.3 has been calculated. In the present example the simple user model is defined by the mean 

Uw = 0m 

and the covariance 

i.e. there is no knowledge available about a preferred direction of the movement of the user. 
The covariance matrix C w only represents the distribution of the step lengths of the user. 
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EKF filter (Extended Kalman Filter) 

[0134] When the EKF filter is used the measurement equation (5) is linearized at the mean of 
the a-priori density function. Fig. 5 shows the results of the EKF for the first and the 170th filter 
step. For each filter step the measurement of the logarithmic receive signal power of a single 
transmitter is used for updating the estimated position. 

[0135] It can be seen that after the first filter step there is a significant estimation error. The 
true position of the receiver does not lie within the support values of the estimated a-posteriori 
density function, which is shown by the profile representation. Furthermore there is no crossing 
point between the true measurement equation and the estimated density. It is thus not possible 
to apply data validation processes to establish whether a measured value can be "explained" by 
the estimated density. 

PDSME (Prior Density Splitting Mixture Estimator) 

[0136] Fig. 6 shows the a-posteriori densities for the first and 170th filter steps estimated by 
the PDSME algorithm. The approximation of the non-Gaussian a-posteriori density in the first 
filter step which lies very close to the numerically calculated optimum Bayesian solution of the 
filter step is clearly much better than in the case of the EKF filter. Data validation procedures 
can thus be successfully applied. 

[0137] In Fig. 7 the expected values of the a-posteriori density functions calculated with the 
EKF and the PDSME are compared to the true position. If this figure is compared with Fig. 8 it 
can be seen that for almost linear filter problems, i.e. if the covariance of the estimated position 
by comparison with the non-linearity of the measurement equation after a number of filter steps 
is small, EKF and PDSME yield almost the same results. In these cases the PDSME uses only 
a modest number of Gaussian mixture components whereas for strong non-linearities at the 
beginning of the localization experiment a higher number of approximation components is 
necessary to reduce the estimation error. The superior power of the PDSME filter by 
comparison with the EKF is also shown by the average estimation error 
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[01 38] over the N F = 1 70 filter steps which amount to 3,30 m for the EKF and 1 ,22 m for the 
PDSME. 

Conclusions 

[0139] In this work a stochastic approach to a localization of radio communications devices 
has been offered to which is based on the measurement of the logarithmic signal power of the 
receivable transmitters by a mobile station. For each transmitter a stochastic measurement 
model has been identified which consists of a deterministic and a stochastic component. This 
measurement model has been used for the estimation of the position of the receiver through an 
innovative Gaussian mixture estimator which is based on the splitting of the a-priori density 
corresponding to a linearization error criterion. This criterion is very similar to the Kullback- 
Leibler distance between the true and the approximated a-posteriori density which is calculated 
by a linearization of the measurement equation. By applying this estimation procedure to the 
localization of DECT mobile telephones significant improvements to the estimation quality can 
be achieved if the PDSME is used instead of the standard approaches such as the EKF. 
Further improvements to the localization quality by identifying better deterministic measurement 
models and more exact identification of the measurement noise can be achieved which do not 
assume any independence between the different uncertainties mentioned in this work. 

[0140] The invention has been described in detail with particular reference to preferred 
embodiments thereof and examples, but it will be understood that variations and modifications 
can be effected within the spirit and scope of the invention covered by the claims which may 
include the phrase "at least one of A, B and C" as an alternative expression that means one or 
more of A, B and C may be used, contrary to the holding in Superguide v. DIRECTV, 
69 USPQ2d 1865 (Fed. Cir. 2004). 
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